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ABSTRACT 

We present new [/-band photometry of the magnetic Helium-strong star a Ori E, obtained over 2004- 
2009 using the SMARTS 0.9-m telescope at Cerro Tololo Inter- American Observatory. When combined 
with historical measurements, these data constrain the evolution of the star's 1"?19 rotation period 
over the past three decades. We are able to rule out a constant period at the p nu n = 0.05% level, and 
instead find that the data are well described (p nu ii = 99.3%) by a period increasing linearly at a rate 
of 77ms per year. This corresponds to a characteristic spin-down time of 1.34Myr, in good agreement 
with theoretical predictions based on magnetohydrodynamical simulations of angular momentum loss 
from magnetic massive stars. We therefore conclude that the observations are consistent with a Ori E 
undergoing rotational braking due to its magnetized line-driven wind. 

Subject headings: stars: chemically peculiar — stars: individual (HD 37479) — stars: magnetic field 
— stars: massive — stars: mass-loss — stars: rotation 



1. INTRODUCTION 

The helium-strong star a Ori E (HD 37479; B2Vpe; 
V = 6.66) has long been known to harbor a circum- 
stellar magnetosphere in which plasma is trapped and 
forced into co-rotation by the star 's strong (~ 10 kG) 
dipol a r magnetic field ( see, e.g., lLandstreet &: Borral 
119781: iGroote k Hungerl I1982D . This magnetosphere 
is largely responsible for the star's distinctive eclipse- 
like dimmings, which occur when plasma clouds tran- 
sit across the stellar disk tw ice every L d 19 rotation cycle 
(Townsend & Owocki 2005). Some fraction of the star's 
photometric variations likely also arise from its photo- 
spheric abundance inho mogeneities, as in othe r chemi- 
cally peculiar stars (e.g. iMikulasek et al.l l2009h ; but for 
a Ori E the mag netospheric con tribution to the varia- 
tions is dominant (|Townsendll2008l ) . 

This paper presents new [/-band photometry of the 
star's primary light minimurrQ, obtained over four sea- 
sons spanning 2004-2009 using the SMARTS 0.9-m 
telescope at Cerro Tololo Inter-American Observatory 
(C TIO). When c ombin ed with historical measurements 
bv iHesser et al.1 (|1977j ). these new data allow a precise 

1 The term 'primary' stems from earl y mis-identifications of the 
star as an eclipsing binary system (e.g., [Hcsscr et al.iri976h : here, 
it simply indicates the deeper of the star's two light minima. 



measurement of the star's rotation period, and its evolu- 
tion, over the past three decades. 

A description of the observations, both archival and 
new, is provided in the following section. In <j3l we dis- 
cuss a procedure for accurately measuring the times i m i n 
of primary light minimum, and then use these measure- 
ments to construct a standard observed-minus-corrected 
(O — C) diagram for the star, allowing us to assess the 
evolution of the star's rotation period. We discuss and 
summarize our findings in |J4] 

2. OBSERVATIONS 
2.1. 1977 

IHesser et alj (|1977f ) observed a primary light minimum 
of a Ori E on the night of 1977 January 26/27, as part 
of their long-term Stromgren uvby monitoring of the star 
using the number 1 0.4-m telescope at CTIO. Their pho- 
tometric data were kindly provided to us in electronic 
form by Prof. C. T. Bolton. No error estimates were sup- 
plied, so a measurement error Am = 7 mma g is assumed 
for all data points — the value quoted by IHesser et al.1 
(1976) as an upper limit on their photometric uncertain- 
ties. We do not corre ct for the apparent s eason-to-season 
brightening noted bv IHesser et al.1 (fl97 7*PI. since this has 

2 This brightening likely stems from long-term variations in the 
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Fig. 1. — Light curves for the primary minima across the five seasons ([01; filled symbols denote those points that contribute to the 
minim um fitting, and open symbols those that do not. Vertical bars indicate the estimated measurement errors, and the best-fit parabolae 
([EHJ arc drawn as the solid curves. 



no effect on the t m i n determinations. The u-band light 
curve is plotted in Fig.Q] the accompanying vby data are 
not shown, since they play no direct role in the period 
determination (however, see 



2.2. 2004 

We observed a primary light minimum on the night 
of 2004 November 26/27, during Johnson UBVRI 
monitoring of a Ori E using the Cassegrain-focus 
Tek 2048 CCD on the S MARTS 0.9-m telescope (see 
lOksala k Townsend! [20071 for the full dataset). For all 
[/-band exposures (in every case, 420s long), the Tek 
#1 filter (center 3575 A, FWHM 600 A) was used in tan- 
dem with the ND3 neutral density filter (7.5 mag at- 
tenuation). CCD frames were reduced using standard 
iraf tasks for bias subtraction and flat-fielding, and 
cosmic rays were cleaned via Laplacian edge detection 

comparison star HR 1861; see Olsen 1 119771 ). 



(|van Dokkum|[200l[ ). The phot task of the DAOPHOT 
package, with an aperture radius of 10 pixels (= 3.96") 
and a sky annulus radius of 15-20 pixels, was used to 
perform synthetic aperture photometry. Measurement 
errors AU were estimated as the sum in quadrature of 
the photometric noise reported by phot, the CCD read 
noise, and the atmospheric sc intillation noise des cribed 
by the expression on p. 141 of IBirnev et al.l (|2006D . 

As a comparison star, we adopt the nearby (30") 
a Ori D (HD 37468D; B2V; V = 6.62) due to its 
color and brightness similarity; the Ub — Ub differential 
light curve is plotted Fig. Q] To make an independent 
check on the constancy of a Ori D, we compare its U- 
band data against a Ori AB (HD 37468; O9V+B0.5V; 
V = 3.80). Over the night, we find a mean differ- 
ence Ub — Uab = 3.240 mag and a standard deviation 
o~{Ub — Uab) = 2.8mmag; the latter is consistent with 
the mean measurement error A{Uu — Uab) — 2.6mmag. 
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2.3. 2006 

We observed a primary light minimum on the night 
of 2006 January 30/31. The one modification to the 
0.9-m instrumental setup was the removal of the ND3 
neutral density filter from the light path, thereby short- 
ening exposure times to 2 s (with the idea of obtaining 
more measurements during the night). In hindsight, this 
was a counterproductive move. Although the proximity 
of a Ori D and a Ori E mean that shutter corrections 
are unimportant, the exposures are strongly affected by 
scintillation noise, especially toward the end of the night 
when the airmass becomes large. Moreover, a Ori AB 
is saturated in all CCD frames and cannot be used as a 
check star. 

The observations were reduced and analyzed as before; 
the Ux) — Ue light curve is plotted in Fig. [TJ 

2.4. 2008 

We observed a primary light minimum on the night 
of 2008 November 25/26. To obtain more-reasonable 
(22 — 65 s) exposure times than in the 2006 season, we 
reintroduced a neutral density filter (ND2; 5.0 mag). 
Changes to the default telescope configuration led to 
the inadvertent substitution of the Tek #2 [/-band fil- 
ter (center 3570 A, FWHM 660 A) in place of the Tek 
#1 filter used previously, but this should have negligible 
impact on our results. Due to a further oversight, the 
CCD gain was adjusted from the previous setting of 3.1 
e~/ADU, to 0.6 e~/ADU; this had the unfortunate effect 
of significantly elevating the readout noise. 

The observations were reduced and analyzed as be- 
fore; the Ue — Ue> light curve is plotted in Fig. [TJ Com- 
parison of a Ori D against a Ori AB reveals a mean 
Uu — Uab = 3.236 mag and a standard deviation o{Ue> — 
Uab) — 12mmag. The latter value is rather larger than 

the mean measurement error A(Uu — Uab) — 9.8mmag, 
but not overly so. 

2.5. 2009 

We observed a primary light minimum on the night of 
2009 November 28/29. We kept the ND2 filter from the 
2008 season, but reverted back to the original Tek #1 [/- 
band filter, and set the CCD gain to 1.5 e~/ADU; this 
led to exposure times between 130 and 250 s. The obser- 
vations were reduced and analyzed as before; the Ue — Ue> 
light curve is plotted in Fig. [TJ Comparison of a Ori D 
against a Ori AB reveals a mean Uu — Uab — 3.224 mag 
and a standard deviation <j(Ue> — Uab) = 4.2 mmag, 
the latter con sistent with the mean measurement error 
A([/ D - Uab) = 4.0 mmag. 

3. ANALYSIS 

3.1. Minimum Measurements 

The determination of the time of light minimum (or 
maximum) in astronomical object s has received signifi- 
cant attention in the literature fse elSterkenll2005|. and ref- 
erence s therein). Historically the [Kwee fc van Woerdenl 
(1956) algorithm had long been the standard tool, but 
more recently polynomial fitting has emerged as a pow- 
erful and robust approach. Thus we determine the time 
of light minimum t m - m for each of the five light curves 
plotted in Fig. [TJ by an adaptive parabola fitting proce- 
dure: 



TABLE 1 

Light-curve minima from parabola fitting 



Season 


Band 


<min (d) 




E 


1977 
2004 
2006 


u 
U 
U 


2443170.6097 t^oon 
2453336.7789 t°mii 
2453766.6753 l^ooog 


1.01 

1.56 
1.11 


329 
8866 
9227 


2008 
2009 


U 
U 


2454796.7624 t^orao 
2455164.7343 l^oo? 


1.12 
0.75 


10092 
10401 



1. The time of the dimmest point in the light curve is 
chosen as the initial i m i n . 

2. Weighted % 2 minimization is used to fit a parabola 
to those points lying within 2 hours of the current 
tmin (this time interval is chosen to exclude the 
non-eclipse parts of the light curve). 

3. A new i m j n is chosen by analytically evaluating the 
minimum of the fitted parabola. 

4. Steps (ii)-(iii) are repeated until the value of i m ; n no 
longer changes. (Typically, 3-4 of these iterations 
are required.) 

The parabolic fits are drawn over the light curves in 
the figure, and Table [TJ documents the i m ; n and reduced 
X 2 values associated with each fit. Here and throughout, 
times are expressed as heliocentric Julian dates (HJD). 
The tabulated t min error bounds a re calculated from 
bootstrap Monte Carlo simulations ()Press "et~aT1[l99l) . 
Specifically, for a given light curve comprising n points, 
a synthetic light curve is generated by selecting n points 
at random with replacement. The parabola fitting pro- 
cedure is then used to determine / m ; n for the synthetic 
curve. This sequence is repeated many times to build up 
a population of synthetic Z m ; n values, from which we de- 
rive a probability distribution function -F^mm) reflecting 
a best estimate of the distribution from which the actual 
tmm measurement is drawn. 

Figure [5] plots the probability distribution functions 
(PDFs) associated with the five i m i n measurements; each 
is based on a population of 10 7 synthetic light curves. 
The error bounds quoted in Table [TJ are the confidence 
limits relative to i m j n that each enclose 34.1% of the area 
under the PDF. Although this choice mirrors the 1-a lim- 
its of a Gaussian distribution, the PDFs in Fig. [5] under- 
score that the analysis here does not presume Gaussian 
errors. Moreover, while the final error in £ m j n is insensi- 
tive to the estimates made for the measurement errors on 
the individual photometric data, the Monte Carlo sim- 
ulations used to derive the PDFs do properly allow the 
actual errors of the noisier data, e.g. in 2008, to produce 
a larger uncertainty in i m ; n . 

3.2. O-C Fitting 

We apply t he standard O — C diagram technique (e.g., 
to assess the evolution of the rotation pe 



we apply t n 
Sterken 2005) 



riod. With the reference epoch defined by 
1976[) and the rotation period measured by 



Hesser et al 



Hesseretal 



19771 ). an observed- minus-calculated value 



o-c = u 



(2442778 d 819 + l d 19801£) (1) 
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Fig. 2. — Probability distribution functions for the time of primary minimum; the thick vertical line indicates the measured t m i n . The 
shaded regions to the left and right of the line each constitute 34.1% of the area under the PDF, and define the lower and upper error 
bounds on t m i n (shown as dotted vertical lines). 



TABLE 2 
Fits to the O — C data 



Fit Type bi x 10 2 (d) b 2 x 10 5 (d) b 3 x 10 9 (d) p null (%) 



Linear 0.38 ±0.12 2.78 0.05 

Quadratic 1.00 to.il 1-29 t°'.io L44 ^o'.ii 99 - 3 



is evaluated for each t m - m in Tab.[TJ here and in the table, 
the cycle number E is the integer that minimizes \0 — C\. 

Figure [3] plots O — C as a function of E; the error bars 
on each point are taken from the corresponding error 
bounds on t m i n . Also shown in this O ~ C diagram are 
fitted linear and quadratic models of form 

(O - C) mod = 61 + b 2 E (2) 



and 

(0-C) mod = b 1 +b 2 E + b 3 E 2 , (3) 

together with the associated fit residuals. The linear 
model corresponds to a constant rotation period, while 
the quadratic one represents a period that increases lin- 
early with time. The coefficients {bj} are determined by 
a non-% 2 maximum likelihood estimation, which seeks to 
minimize the statistic 

Q = J2-^gF[o-c-(o-CT od ]; (4) 

here, the summation is over the five light minima, each 
with its respective O-C, (O - C) mod and F{t min ). For 
a given model, the quantity exp(— Q) is proportional to 
the likelihood that the measured O — C values could have 
arisen by chance-fluctuation departure s from the model . 
We use a downhill simplex algorithm ([Press "et~all[l992l) 
to minimize Q, iterated until the maximum relative de- 
viation between simplex vertices drops below 1CP 4 . 
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Fig. 3. — The observed-minus-corrected diagram for the primary 
minimum measurements. The solid (dashed) lines indicate the 
best-fit quadratic (linear) models; the residuals relative to these 
models are shown below the O — C diagram as filled (open) sym- 
bols. 

Table [2] summarizes the coefficients {bj} of the linear 
and quadratic models. As in £13.11 error bounds are de- 
termined via Monte Carlo simulations. However, rather 
than generating synthetic O — C data by bootstrapping, 
we construct them by perturbing each i m j n with random 
deviates drawn from the appropriate PDF (cf. Fig. [5]). 
From the same simulations we also determine the dis- 
tribution of the Q statistic, enabling us to associate the 
linear and quadratic fits in Fig. [3] with a likelihood p n uii 
(also specified in Tab. [2J that the null hypothesis is true: 
that is, that the deviation of the O — C values from the 
model arises purely due to chance fluctuations. 

3.3. Period Evolution 

Table [2] indicates that the null hypothesis for the linear 
O — C model is extremely unlikely (p nu ii = 0.05%), al- 
lowing a constant rotation period to be ruled out with a 
high degree of confidence. However, the converse is true 
for the quadratic model, which fits the data extremely 
well. Combining its coefficients with eqns. (fTJ) and @ 
gives a revised ephemeris for the primary light minimum 
of a Ori E as 

t min = 2442778 d 8290i° : ^l + 

l d 1908229±°;™^ + ^^-Vii * W- 9 E 2 . (5) 

By taking the derivative with respect to cycle number, 
the instantaneous period is found as 

P = l d 1908229t°«§ + 2 d 89t° : ^ x 1Q-°E, (6) 

which grows linearly at a rate P — 77 ms per year. (We 
discuss the implicit assumption of smooth period growth 
in []4]). At the reference epoch E = 0, this period is 
ra ther larger than the P = l d 19081 ± d 00001 reported 
bv lHesser et al.l ()1977f j: howeve r these author s ' erro r es- 
timate seems overly optimistic. iReiners et al.l (|2000f ) de- 
termined a period P = l d 19084 ± d 00001 by combining 
new helium-line data with historical measurements from 



iPedersen fc Thomsenl (|1977l ): their value is in good agree- 
ment with the mean period P = l d 190833 obtained by 
averaging eqn. © over the 1976-1998 interval spanned 
by the helium data. 

3.4. Systematics 

Before discussing the significance of the measured pe- 
riod increase, we briefly review factors that may have 
a systema ti c effe ct on this result. As demonstrated by 
iTownsendl (<2008[ ). the timing of light minima is sensi- 
tive to the optical depth of the magnetospheric plasma 
clouds. In principle, the progressive shift in t m i n toward 
later times (seen in the residuals plot of Fig. [3]) could 
be explained by the secular accumulation of plasma i n 
the magnetosphere (see, e.g. JTownsend fc Owockil2 005) . 
Howev er, based on the examples given by ITownsendl 
(2008, his section 2.4), the observed shift of « d 02 be- 
tween the 2004 and 2009 seasons would require a factor- 
six increase in the optical depth of the magnetosphere. 
The minima in Fig. [T] clearly do not exhibit this kind 
of dramatic change. Indeed, although there are some 
season-to-season variations in the minima depths (on the 
order of 20 mmag) , no monotonia trend is seen; we be- 
lieve the variations are probably due to changes in the 
distribution of scattered light from a Ori AB. Accord- 
ingly, we rule out the possibility that the i m i n shift is 
due to plasma accumulation. 

Similar reasoning can be addressed to concerns over 
the change in filters, from Stromgren u in 1977 to John- 
son U in the later observations. In passbands where the 
magnetosphere is more opaque, the time of primary light 
minimum will tend to occur later. This color-fm i n cor - 
relation can be clearly seen in the iHesser et al.l ()1977f ) 
observations; in the Stromgren u band (falling blueward 
of the Balmer jump), the time of primary light minimum 
is d 007 later than in the v band (falling redward of the 
Balmer jump). Because the U band straddles the Balmer 
jump, the expected time lag between u and U should be 
around half of this, ~ d 003. An adjustment of this or- 
der to the 1977 t m i n point, to correct for the color-t m ; n 
correlation, has a negligible effect on our results. 

A final possible issue comes from the use of parabolae 
to measure the p rimary minimum times. As discussed 
by iSterkenl (|2005f ). quadratic fitting is often eschewed in 
light-curve analysis on the grounds that it is unable to 
adequately model asymmetric light minima. However, 
if the shape of the light curve does not vary from cycle 
to cycle, then this bias is irrelevant: it matters not that 
^min occurs slightly before or after the precise time of 
minimum, as long as the lead or lag remains invariant. 

Nevertheless, to explore any bias introduced by the use 
of quadratic fitting, we have repeated our analysis using 
cubic fitting to measure the i m i n values. Three salient 
points stand out from this re-analysis: (i) the x^ ed values 
of the cubic fits are not significantly smaller than those 
of the quadratic fits (cf. Tab. [I]); (ii) there's no evidence 
for a systematic lag or lead between the quadratic and 
cubic minima; (iii) there is an obvious difference between 
the widths of the PDFs, which are a factor ~ 2 broader 
in the cubic cases than in the quadratic ones. Thus, 
we conclude that quadratic minimum fitting does not 
introduce any appreciable bias, and moreover is the more 
robust approach. 
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4. DISCUSSION 

The Helium-s t rong star HD 37776 was found by 
Mikulasek et al. (2008) to exhibit a progressive length- 
ening in its r?5387 rotation period, with a character- 
istic spin-down time r sp i n = P/P = 0.25 Myr. For 
a Ori E the absence of photometric data in the 1980s and 
1990s means that we cannot empirically differentiate be- 
tween steady spin-down and a sequence of abrupt braking 
episodes. However, the steady scenario is lent strong sup- 
port by magnetohydrodynamical (MHD) simulations of 
angular mom entum loss in magneti cally channelled line- 
driven winds (jud-Doula et al.ir2009( ) , which indicate that 
the lengthening of rotation periods should be a smooth 
process. Therefore, the use of a quadratic ephemeris 
(cf. eqn. [5]) appears justified, and we derive a charac- 
teristic spin-down time T sp i n = 1.341q 'Jg Myr. This value 
coincides very well with t he T sn in = 1.4 Myr pre dicted 
specifically for a Ori E by ud-Doula et al.1 p00 9). from 
their MHD-calibrated scaling law for spin-down times. 
(Such a close agreement is partly fortuitous, given the 
uncertainties in stellar and wind parameters). Assuming 
that T S pin has remained constant over the lifetime of the 
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